Luminosity, selfgravitation and nonuniqueness of stationary accretion 
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We investigate newtonian description of accreting compact bodies witii liard surfaces, including luminosity 
and selfgravitation of polytropic perfect fluids. This nonlinear integro-difi'erential problem reduces, under ap- 
propriate boundary conditions, to an algebraic relation between luminosity and the gas abundance in stationary 
spherically symmetric flows. There exist, for a given luminosity, asymptotic mass and the asymptotic tempera- 
ture, two sub-critical solutions that bifurcate from a critical point. They difi"er by the fluid content and the mass 
of the compact centre. 



INTRODUCTION. 

We study in this paper a newtonian accreting systems con- 
sisting of a compact body enclosed by a ball of selfgravitating 
gas. Our aim is to investigate the following inverse problem: 
assume that one knows the luminosity, total mass, asymptotic 
temperature (or, equivalently, the asymptotic speed of sound) 
and the equation of state of the gas. Assume that the gravita- 
tional potential has a fixed value at the boundary of a compact 
body, but the radius of the body is arbitrary to some extent. 
The question arises; can one specify the mass of the compact 
body? And if one can, is the specification unique? That issue 
is interesting in itself within the framework of astrophysics 
but it is also related to following fundamental physical aspect. 

The standard point of view that compact objects of a mass 
exceeding several solar masses constitute black holes has been 
challenged by Mazur and Mottola 1 1], who developed the no- 
tion of gravatars — stars almost as compact as black holes but 
dispossessed of event horizons. Gravastars are supposed to 
be built from matter which violates energy conditions [2] and 
therefore they can avoid the standard compactness limits (like 
the Buchdahl theorem 831]) or mass bounds. The concept of 
gravastars originated from the belief that the only place after 
Big Bang that quantum gravity could manifest is the collapse 
of matter to a pre-(event horizon) phase, and that quantum 
gravity effect may or can actually stop the collapse. Yet an- 
other conclusion, if one believes in this argument, is that the 
search for quantum gravity should start from the revision of 
the scenario of the gravitational collapse to a black hole. 

The issue of the validity of gravastars to astrophysics can be 
verified, at least in principle, by observations. The fundamen- 
tal question is whether observations can distinguish compact 
bodies endowed with a hard surface (gravastars) from black 
holes. Narayan presented observational arguments in 
favour of the existence of black holes. Abramowicz, Kluzniak 
and Lasota [5] raised several objections, expressing their view 
that present accretion disk models are not capable to make 
distinction between accretion onto a black hole or a gravastar. 

We choose below the simplest selfcontained model that can 
be interpreted as a radiating system. While we do not address 
directly the problem of distinguishing between a gravastar and 
a black hole, our results show that there is an ambiguity even 
in the newtonian level. There can exists at least two solu- 
tions (two compact stars with a hard surface) to the inverse 



problem described above. This is not quite unexpected. In 
a recent general-relativistic analysis (valid also in the newto- 
nian limit) the mass accretion rate M behaves like ^^(1 - y), 
where y = M^/M is the ratio of (roughly) the mass of the 
compact core to the total mass [6]. Thus there exists a maxi- 
mum of M when y - 2/3. The luminosity equals to the mass 
accretion rate M multiplied by the available energy per unit 
mass — the potential 0(/?o) at the hard surface of a compact 
star. Hence this simple analysis of accretion would suggest 
the existence of two weakly luminous regimes: one rich in 
fluid with nif x M and the other with a small amount of 
fluid, nif/M <g: 1. To a given luminosity x might have cor- 
responded two systems with different y's. The actual situation 
in the case with the luminosity is less clear, because the lu- 
minosity impacts also the accretion rate M and there exists a 
complex algebraic relation y = y(x), but the nouniqueness is 
still retained. 

THE SHAKURA MODEL. 

Stationary accretion of spherically symmetric fluids with 
luminosities close to the Eddington limit has been investigated 
since the pioneering work of Shakura |8]. In |8] the gas pres- 
sure and its selfgravity have been ignored and the analysis has 
been purely newtonian. Later researchers included both these 
aspects |9] and extended the analysis onto relativistic systems 
( filOJ — [121). We will consider a spherically symmetric com- 
pact ball of a fluid falling onto a non-rotating compact body. 
The areal velocity f/ of a comoving particle labelled by coor- 
dinates (r, f) is given by U (r, f) - diR, where f is a comoving 
(Lagrangian) time, p denotes pressure, L(R) and Le are the lo- 
cal and Eddington luminosities, respectively. The quasilocal 
mass m{R) is defined by dRm(R) = AnR^g. The mass accretion 
rate is 

M = -AnR^gU. (1) 

g is the baryonic mass density. The equation of state is 
taken to be polytropic, p = Kg^ with F being a constant, 
(1 < r < 5/3). We assume a steady collapse of the fluid, 
which means that all its characteristics are constant at a fixed 
areal radius R: c?/X|s=const = (di - {dtR)dR)X = 0, where 
X - g,U,a^ . a - -\Jdpp is the speed of sound. Strictly saying, 
a stationary accretion must lead to the increase of the central 
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mass. This in turn means that the notion of the "steady ac- 
cretion" is approximate — it demands that the mass accretion 
rate is small and the time scale is short, so that the quasilo- 
cal mass m{R) does not change significantly. The total mass 
m{Rao) will be denoted by M. Furthermore, we assume that 
the radius R^o of the ball of fluid and other boundary data are 
such that Ul « '-^ « a^. 

The steadily accreting gas is described by a system of 
integro-(ordinary)-differential equations. They consist of the 
(Euler) momentum conservation equation 



UOrU = - 
the mass conservation 



Gm(R) 
'~R^ 



-ORP + a—, 

Q R^ 



and the energy conservation 



Lo - L{r) = M 



r- 1 



fl2 f/2 
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The constant a is defined by a = cr/ {AnnipC^, Lq is the total 
luminosity and <p(R) is the newtonian gravitational potential. 
Strictly saying, in Eq. dUl should appear ^ 
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this is well approximated by due to the aforementioned 
boundary conditions. 

We define the Eddington luminosity as that at the outermost 
layer of the accreting gas, that is as = GM/a, while the to- 
tal luminosity Lq is equal to the product of the mass accretion 
rate and the total available energy per unit mass 0o s |(^(/?o)|, 
where Rq is an areal radius of the boundary of the compact 
body, Lo - M<po ([8]). The quantity <po is fixed and Ro in 
the constructed configurations is by definition the areal radius 
at which the absolute value of the surface potential equals to 
(po- Dififerentiating M with respect to R, using the equality 
dRa^/(r-l) = a^dRlng - Orp/q and combining the obtained 
equation with ^ yields the diff'erential equation c?« In L - ^ . 
Its solution can be written in terms of the Eddington and total 
luminosities as follows 



L - Lq exp 



-LqRq 
Lf.R 



(5) 



Here = GM l\(f>{RQ)\ is a kind of modified size measure of 
the compact body. In the case of test fluids = Ro. In the 
case of gravastars one should assume that \(p(Ro)\ ~ and 
^0 ~ IGMjc^ is the Schwarzschild radius; while gravastars 
would require the general-relativistic treatment, we believe 
that essential features will be the same as in the newtonian 
description. 



have been obtained for subsonic solutions. The sonic hori- 
zon (sonic point) is a location where \U\ - a. In the following 
we will denote by the asterisk all values referring to the sonic 
points, e.g. «», f/», etc. Differentiation with respect the areal 
radius will be denoted as prime '. One finds from the mass 
conservation U' - -U{(g' /g) + 2/R). Inserting that into the 
Euler equation (|2]l, one finds 



g R\ 2R 



Gm(R) L(r)a\ 



2R 



(6) 



thus at the sonic point the three characteristics, a,, U, and 
M»/R» are related. 



Gm* I L,a 

al = ul = -TT^ 1 1 - 



2R^ \ GM, 



Gm* I L^M 



LeM, 



(7) 



It is clear from (jT) that in order to have a critical flow one 
should have 



LeM, 



< 1. 



(8) 



It is useful to express the velocity U and the mass density g in 
the following way 



Rllal^"''-'' 



where f/, is the negative square root, and 



(9) 



(10) 



where the constant g^o is equal to the mass density of col- 
lapsing fluid at the boundary Roc- There exist four stationary 
branches crossing at the sonic point, one of which is the ac- 
cretion flow. 
Define 



Lo 



Ro 



^ Le' ^ M' ^ R,' 
K=Lo-U- 2Mal ' 



yLE - L, 



GM, 
Rt 



(11) 



In new variables the necessary condition ([8]l for the existence 
of a sonic point reads .Ji:exp(-X7) < y. A straightforward 
algebra leads to the following equation 



2 \M 



(12) 



THE MASS ACCRETION RATE. 



We shall assume that y < 1 and xy ^ \. From Eq. O follows 
now 



In what follows it is assumed that there exists a sonic point, 
but we believe that the same (nonuniqeness) result would 



Ut = — - 1 - - exp(-xr) 
2R, \ y 



(13) 



3 



Then A* = up to the first order while *P, = AGn drrg. 
One can show, using arguments similar to those applied later, 
that «: al. Therefore (fT2l) simpUfies to 



(14) 



5 -3r 



that is exactly the same result as in the Bondi model [,13]. A 
careful reader would notice that the above derivation requires 
that r is somewhat smaller than 5/3. 

The rate of the mass accretion M as given by ([T]) can be 
expressed in terms of the characteristics of the sonic point. 
One obtains, using equations (|7]l and ([T]i, that 



^2 \ 2(r-i) 



(15) 



Now, one can can show in a way similar to that applied in the 
case of relativistic accretion (^) that under the previously as- 
sumed conditions and >> one has = Xca (M - M,) - 
Mxoo (1 - y) for r 6 (1, 5/3 - 5) with some small 5. 

More specifically, one needs to show that the energy con- 
servation equation (|4|l yields 



aHR) 



<i-(r-i) 



(16) 



two solutions y(x,/3) for any parameter yS, < y < 1 and the 
relative luminosity x smaller than a certain critical value a. We 
shall add the adjective "critical" to all characteristics of this 
critical point; hence a critical luminosity is that corresponding 
to X = a. The obtained results are following. 

Theorem. Define F{x,y) = x-/3(l -y)(y - xexpi-xy))^. 
Then 

i) there exists a critical point x - a,y = b of F, that is 
F{a,b) - and dyF(x,y)\a,b — 0, with a and b satisfying 

0<fl</7<land/7 = Z+'-expC-.y) ^^ ^ ^^^^ _ ^^3. 

ii) For any < x < a there exist two solutions y{x)t bifur- 
cating from {a,b). They are locally approximated by formulae 



y 



J(a - x)(b + a exp(-fl'y)(l - 2ay)) 
+ -/7+- ^ ^ ^— . (20) 

-iJPib - a exp(-fly))(l - a exp(-ay)) 



iii) The relative luminosity x is extremized at the critical 
point (fl, b). 
Proof. 

The two criticality conditions yield two equations 



a-/3(l- b)(b -aexp(-ay)f = 0, 
b - a exp(-fl7) = 2(1 - b). 



(21) 



From Eqs (1211 1 one immediately obtains the desired expres- 
sions 



bounding \(p(R)\ from above by GM/R one obtains 

^<i.(r-i).^^ 



Rai 



(17) 



\i/(r-i) 



Therefore ( fTOl i implies g < Qoa (^l + ^f) • The conclu- 
sion follows now from using this bound in M - M, = dVg 
and appealing to the asymptotic conditions laid down in the 
forthcoming text. The proportionality constant is roughly 
the inverse of the volume of the gas outside of the sonic 
sphere. Since the total luminosity Lq = ip{Ro)M and al/a^ 
is given by ( fT4b . one obtains that 



Lq = (poG TT 



(1 -y){y - xexp(-xy)) 



(5-3r) 
2 \ 2(r-i) 



5-3r 



(18) 



THE ANALYSIS OF THE LUMINOSITY EQUATION. 



2H-aexp(-fly) , 
b = ^- a = 4yS(l - bf. (22) 



Inserting the second of eqs ( l22l i into the second equation in 
(ISTT i. one arrives at 



^ = ^ + y (1 - bf exp(-4/3(l - bfy). 



(23) 



Both sides of this equation are continuous functions of b and 
at = the right hand side of ( |23] | is bigger than the left 
hand side, while alb - I the opposite holds true. Therefore 
there exists a solution. A closer inspection shows that there 
is only one critical solution. This fact actually guarantees that 
solutions y(x) that bifurcate from {a, b) extend onto the whole 
interval x e (0, a). (Hint: use the implicit function theorem 
and the fact that (a, b) is the unique critical point). The specific 
form of a solution close to a critical point can be obtained as 
follows. Insert x - a + e,y - b + y{e) into F{x,y) - and 
expand F keeping the terms of lowest order One obtains a 
reduced algebraic equation 



It is convenient to cast the formula (fTST i in terms of the rel- 
ative luminosity, x - LqILe, 



x^l3{\-y){y-xexp{-xy)) , 



(19) 



(5-30 

2(r-i} 



where l3^X^4>^aGn^^ [^Y' 

We shall analyze solutions of this cubic (with respect to the 
variable y) equation. One can show that there exist at least 



{\+2l5{\-b){b- ae-"^) {ae-^^y - e-"^)) e + 
j3(3b-l-2ae-''^)y^(e) = 0. 



(24) 



Notice that at the critical point 3b - I - 2a exp(-fly) = 1 - 
a exp(-fly) (see the second of Eqs (|2TI) ) while 

2jS (1 - b){b - a exp(-a7)) {a exp{-ay)y - exp(-fly)) - 
2a {a exp{-ay)y - exp(-fly)) jib - a exp(-fly). (25) 
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Inserting (IZST i into ( |24] | and finding ^(e) from the latter imme- 
diately leads to the approximate solution of ii). 

Now we shall prove the third part of the Theorem. Let 
(xo,yo) be a non-critical solution of F{x,y) - with the do- 
main being subset xe^" < yof the square 0<jic< 1,0 <y < 1. 
Then from the implicit function theorem there exists a curve 
x(y) such that F(x(y),y) - for y belonging to some vicinity 
of yo. Along this curve one has 

dx _ (y- xe-^y) (Sy - 2 - xe-'^) 

Ty "^l+2/?(l-y)(y- 



Bifurcation point =0.570, 0.8258616873 



■ (26) 

xe-^y^e-^y - jxe-'^y) 



The denominator is strictly positive while the nominator van- 
ishes only at critical points. That proves the assertion. It fol- 
lows from the form of the approximated solution constructed 
in part ii) that x = a is the extremal value of the relative lumi- 
nosity. That ends the proof of the Theorem. 

It is easy to observe that the coordinates a, h of the critical 
point increase with the increase of the parameter Indeed, 
from Eq. (i22l i follows ^ = e"°''(l - ay); this is bigger than 
zero, since ay < I. Differentiation of the first equation in (i22l l 
with respect to yields § = 4(1 -bfl{% + 12y6(l -bf) > 0; 



we exploited here the fact that ^ > 0. Since along the critical 
curve ^ = ^ obtain also ^ > 0. In the particular case 

of y = one can explicitly find b by solving (i23T l. 



1 



1/3 



(27) 



(i22l i gives then a(J3) and one can check explicitly that both 
a and b monotonically increase with the increase of /3. The 
forthcoming figure shows the two branches xi(y), X2iy) bifur- 
cating from a critical point, in the case when y = 0. 

The next interesting fact is that at critical points the param- 
eter b is not smaller than 2/3. This lower bound is saturated at 
small j6, i. e., when the relative luminosity a is small (notice 
that a < p), which corresponds to the maximal rate of the mass 
accretion in irradiating systems @]. Since the fluid abundance 
is equal io \ - b, we can conclude that critical configurations 
have less fluid than 1/3 of the total mass, and the upper bound 
1/3 is saturated in the limit of vanishing radiation. The max- 
imal mass accretion rate does not occur exactly si y - 2/3, 
as in systems with no radiation |i6J, but at a somewhat larger 
value. 



DISCUSSION. 

In the foregoing text we have assumed the existence of 
an accreting system satisfying particular boundary condi- 
tions. Although the detailed behaviour is given by the integro- 
differential nonlinear system (|2l|4]l, its essential features can be 
deduced from the algebraic equation ( fT9] l. The validity of this 
projection of the original equations to the algebraic problem 
relies on the existence of appropriate solutions. We checked 
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FIG. 1: Here p = 50. The mass abundance is put on the ordinate 
while the relative luminosity x is shown on the abscissa. The circle 
encloses the bifurcation point (a, b). 



numerically that there do exist such solutions and that the rel- 
ative error made in the above approximations is of the order 
of 10"^. Here follows an example. 

The mass of the system is 100 solar masses {M - 1.989 x 
lO-*^ grams), its size 1.254 x lO''* cm and the relative lumi- 
nosity X = !^ - 0.1. The adiabatic index was taken F = 4/3, 
the asymptotic speed of sound = 1.875 x 10^ cm/s and 
,Ro = 6 X 10^ cm. Let us remark that if M and Rq are fixed, 
then also the potential at the surface of the compact core is 
fixed. The set of quantities hitherto specified should in prin- 
ciple determine a solution; any additional information would 
lead to a contradiction. 

We have found, by integrating the original integro- 
differential equations (|2] — |4]i that there exist two solutions 
with following characteristics: i) solution I — /?» = 1.2 x 10^ 
cm, = 0.1 IM and goo ~ 2.144 x 10"'^ g/cm^; ii) solution II 
— = 1.068x10^ cm, M, = 0.9999M and p^o ~2.41xl0"'2 
g/cm^. The only diflFerence between numerical solutions and 
algebraic approximations shows in q^o and it is of the order of 
0.1% . 

Under a simplifying assumption that y <«c 1 one can ne- 
glect the exponential term in xexp(-xy). The analysis of (i22l i 
and (IZTT i allows us to conclude that a can be as close as one 
wishes to 1 for large fi. Putting that in physical terms, the lu- 
minosity L() can be close to the Eddington luminosity Le if the 
parametr y6 is large enough, or more precisely, if the product 
{MIM()f- lifloolaooof (where Mo and aoco are some reference 
quantities) is large enough. The two bifurcation solutions have 
a luminosity x < a and they are characterized by two differ- 
ent numbers y\,y2; the latter can be markedly different (say 
y\ > IQyi) only for x a. This is suggested by the point ii) 
of the Theorem which shows the branching of the two solu- 
tions from a critical point, but numerical calculations confirm 
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that more precisely. A practical conclusion is a follows: if a 
given system radiates with the luminosity close to its critical 
luminosity (and, in particular, to the Eddington limit), then the 
cores corresponding to the bifurcating solutions have similar 
masses. The necessary and sufficient conditions for having 
two cores with significantly differing masses. Mi « M2, is 
that the luminosity is much smaller than Le or the critical lu- 
minosity fl, respectively. 



FINAL REMARKS. 

In conclusion, we have proven that there exist two radi- 
ating systems having the same mass, luminosity, asymptotic 
temperature and surface potential. In extreme cases — for in- 
stance when the total mass is very large and the luminosity low 
— one of the solutions corresponds to a compact body having 
a mass close to the total mass and a small amount of gas, while 
the other solution consists of a light compact body enclosed 
by a heavy cloud of gas. Returning to the question touched at 
the beginning - one would not make a distinction between a 
gravastar and a neutron star within the simple model discussed 
here. 

There is a number of obvious questions that can be asked in 
connection with presented results. 

The first question is whether one can justify the assump- 
tion — which we make — of stationary accretion in the case 
of a system where the cloud of accreting gas is heavier than 
the compact core. The resolution of this problem would re- 
quire the investigation of the dynamical version of the Shakura 
model. The stationary solution of equations (|2] — HI should 
be inserted as initial data to the dynamical equations and if the 
dynamical solutions remains close to the stationary data then 
the assumption would be said to be justified. We performed an 
analogous verification of the assumption of stationarity in the 
case of selfgravitating gas (in the framework of general rel- 
ativity) {(^, f^) with the positive conclusion. The evolving 
system remained essential unchanged in its interior for times 
much smaller than (Rca/aoa)- We believe that the same con- 
clusion will be valid for the Shakura model. 

The second question is — accepting that the approximation 
of steady accretion is correct — whether one can reduce the 
integro-diffferential problem to an algebraic one. Our paper 
provides a proof that this is so if a solution satisfies particu- 
lar boundary data. Numerical examples suggest that the set of 



solutions that satisfy the required boundary conditions is not 
empty. On the other hand, we believe that this nonuniqueness 
of solutions will manifest in a much larger set of steadily ac- 
creting systems which do not necessarily satisfy our boundary 
conditions; that it can be generic. The restriction to a particu- 
lar sample of systems was made only in order to simplify the 
analysis. 

Finally, there arises the question that inspired the paper: 
whether one can unambiguously distinguish between accre- 
tion onto compact bodies from accretion onto a black hole. 
This would require the extension of the above into general- 
relativistic systems. That is, from the general-relativistic point 
view, a straightforward exercise (see for instance 
is, however, a fundamental difficulty ( ifioll — 



14]). There 
1 11) in de- 



scribing interaction between accreting perfect fluid and the 
nonisotropic (and thus nonperfect) outflowing radiation. This 
introduces another level of ambiguity into the problem, but 
we believe that the nonuniqueness, that was discovered in the 
newtonian accretion will remain. 

This paper has been partially supported by the MNII grant 
1P03B 01229. 
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